Thermal Denaturation of Fluctuating DNA Driven by Bending Entropy 
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A statistical model of homopolymer DNA, coupling internal base pair states (unbroken or broken) 
and external thermal chain fluctuations, is exactly solved using transfer kernel techniques. The 
dependence on temperature and DNA length of the fraction of denaturation bubbles and their 
correlation length is deduced. The thermal denaturation transition emerges naturally when the 
chain fluctuations are integrated out and is driven by the difference in bending (entropy dominated) 
free energy between broken and unbroken segments. Conformational properties of DNA, such as 
persistence length and mean-square-radius, are also explicitly calculated, leading, e.g., to a coherent 
explanation for the experimentally observed thermal viscosity transition. 

PACS numbers: 87.10.+e, 87.15.Ya, 82.39.Pj 

Double-stranded DNA (dsDNA) is made up of two 
intertwined interacting semi-flexible single-strand DNA 
(ssDNA) chains. Over fifty years ago it was recognized 
that the intracellular unwinding of DNA at physiolog- 
ical temperature has as counterpart the thermally in- 
duced denaturation above physiological temperature of 
purified DNA solutions where dsDNA completely sep- 
arates into two ssDNA. Despite the differences between 
the two mechanisms, this observation has led to an inten- 
sive study of thermal denaturation [J, Q ■ The stability 
of dsDNA at physiological temperature is due to the self- 
assembly of neighboring base pairs within a same strand 
via base-stacking interactions and of both strands via hy- 
drogen bonds between complementary bases. The bond- 
ing energy is, however, on the order of ksT (thermal en- 
ergy) 01 and thermal fluctuations can lead, even at phys- 
iological temperature, to local and transitory unzipping 
of dsDNA The cooperative opening of consecutive 
base pairs leads to denaturation bubbles and the melt- 
ing temperature, T m , above which bubbles proliferate, 
depends on sequence, chain length, and ionic strength. 
Experiments show, for example, that there exist a bub- 
ble initiation barrier of ~ lOfcsT and free energy cost 
of ~ 0.1/cbT for breaking an additional base pair in an 
existing A-T bubble A detailed understanding of 
equilibrium [l[ and dynamical [B[ properties of DNA in 
solution is still being sought and a consensus concerning 
the physical mechanism behind the denaturation transi- 
tion has not yet been reached. 

A variety of mesoscopic models have been proposed to 
account for the thermodynamical properties of denatura- 
tion bubbles in DNA. They range from i) simple effective 
Ising-like two-state models [H to more detailed ones such 
as ii) loop entropy models (with or without chain self- 
avoidance) 12, la, ll| , and hi) non- linear phonon models, 
where the shape of the interaction potential between base 
pairs is more precisely taken into account @, 0| . To get a 
transition in models i) and ii), an effective temperature 
dependent base-pair chemical potential must be inserted 
by hand. For finite chains, type (ii) models simply refine 



the sharpness of the transition [l[ , but do not attempt to 
provide a deeper explanation of the physical mechanism 
- our aim here. For type (iii) models it has been shown 
that there can be a denaturation transition analogous to 
interface unbinding, due to a gain in translational en- 
tropy. If, however, physically reasonable values for the 
model parameters are used [1, 0, , T m appears to be 
much too high and the transition width much too large. 

It has been shown experimentally that dsDNA is two 
orders of magnitude stiffer than ssDNA at normal salt 
concentration. We show in this Letter that taking into 
account this difference in bending rigidity provides a 
novel physical explanation for the denaturation mecha- 
nism and leads to realistic values for transition temper- 
atures and widths. Despite important recent advances 
in understanding the crucial role played by DNA bending 
rigidity in explaining force-extension [ill Il2| and cycliza- 
tion [13j experiments, its importance has not yet been 
clearly elucidated in the context of denaturation. We 
show, via a well-defined coupled Ising-Heisenberg statis- 
tical model, that an entropy driven denaturation transi- 
tion emerges by integrating out chain fluctuations, due to 
the entropic lowering of the energetic barrier for bubble 
nucleation. This minimal model neglects all other (resid- 
ual) interactions between bases arising from, e.g., elec- 
trostatics [HI, self-avoidance [f|, and helical twist (lit ]. 

We begin by considering a worm like chain (WLC) 
Hamiltonian H[vi,Y2\ for two interacting ssDNA ho- 
mopolymers of length L: 



H 



ds 



|«6(p(«))r?(«) + |«.(p(*))*?(«) 



dsV(p(s))), 



(1) 



where s is the curvilinear index, r^(s) is the position of 
chain i at base position s, and p(s) = ri(s) — r^s) is the 
relative (internal) coordinate (f3 = 1/ksT and r(s) = 
^j). The coefficient Kb is a bending elastic modulus that 
is proportional to the short distance cut-off £q as 0.34 nm 



2 



(the monomer length) and K s (p) = l/[{3 2 n b (p)]. In order 
to account for the enhanced stiffness of dsDNA, K b must 
depend on p(s), e.g., n b (p) = K b s + AK b e~( p ~ p °)/ A , where 
AK b = nf s /2 - Kf > and A ~ lA is the range 0. 
The persistence length of ssDNA is i s p s — (3nl s 1 nm 



and that of dsDNA is l d p s 



PkI 8 w 50 nm (at 300 K 
and physiological ionic strength). The potential V(p) ac- 
counts for both the effective short range hydrogen bond- 
ing interaction between complementary bases at the same 
s and part of the stacking interaction between neigh- 
boring bases; a convenient form is the Morse potential: 
V{p) = (D/£ )v((p - po)/X) where v(x) = e~ 2x - 2e- x , 
leading to a well depth of D/£q at po [8j|. 

At high temperature, p(s) ~ po ^ A, and there- 
fore the system decouples into two semi-flexible non- 
interacting ssDNA chains. After introducing the center- 
of-mass (external) coordinate, X(s) = [ri(s) + r 2 (s)]/2, 
the partition function can be rewritten as a sum 
over center-of-mass and relative configurations: Q = 
f DpZ>Xexp{-/?(ir oxt [X,p] + H mt [p])}, where 



H h 



H cx t — 



. f d S [ Kb (p)p 2 + Ks (p)p 2 + V(p)] (2) 
Jo 



ds 



K b (p)± 2 + K S { P )X 2 



(3) 



We integrate over X to obtain an effective 
model for p, Q = J 2?pexp{— f3Heg [p]}, where 
H cS [p] = H int [p] + F c *t[p] with J"ext[p] = 
— fc^rin [J 2?Xexp {— (3H cxt [X-, p]}] . The external 
free energy ^cxtlp] at frozen p can be evaluated by 
introducing the local center-of-mass tangent vectors 
t = X and then changing variables to t = t/ \/ K b (p). 
There are two contributions to ^"cxt: one will renormalize 
the second term in i?i nt and the other will renormalize 
the potential to Vr. The latter mainly gives rise to a 
purely entropic barrier favoring bubbles, 



VR{p) = V{p) + \ k -^\n 



(4) 



which lowers the well depth from D to Dr w D — 
(3/2)/c B Tln [ref /(2«g a )]. The entropic bending contri- 
bution can be extremely important when D is in the com- 
monly accepted range off < D / 'IzrT < 5 at T — 350 K 
and nf s /(2k^ s ) ~ 25. Scaling arguments then show that 
the melting temperature gets reduced by a factor ~ 2 
down to experimental values 



17] 



In order to illustrate the above mechanism in more 
detail, we introduce a discretized, exactly soluble, version 
of the above model, which captures the essential physics. 
We map the external tangent vector, t(s), to tj (s = 
and an internal variable I — 2 0(p(s) — po — A) (with 
the step function) to an Ising variable: <jj = 1 for 
an unbroken bond (state A) and <n = — 1 for a broken 
one (B). Each link vector can be denoted by the solid 



angle fi. 



i,4>i). The energy H[<Ji, tj] of a state is 



N-l 



H = ^+i^-U-t l+1 )+H l (J,K,p), (5) 
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mE °"< 



i=l 



Thus contains only the dimensionless parameters 
Ki,i+i = J = /3J, p = (5p, and if = f3K. This 

type of model was introduced in the context of helix-coil 
transitions in 2D [3] and later used in various forms to 
study DNA force-extension and cyclization 11 , 13, EH • 
The first term in H, corresponding to H cxt in the con- 
tinuum model, is the bending energy of a discrete WLC 
with a local rigidity Kj^+i = ka = (3nf s /do for a nearest 
neighbor link of type A-A, kr = 2/3k^ s /£o for B-B and 
kar for A-B. In the Ising part, Hi, the first term mim- 
ics the gradient terms in Eq. @ and accounts for the 
local destacking energy Q. The second term accounts 
for the difference in stacking energy between a segment 
of dsDNA and a denaturation bubble. The third term 
corresponds to the energy (2/2) required to break a base 
pair, contributing to D in the continuum model. There is 
evidence that K <C p, which justifies the choice of K = 
adopted below jlij . The bare parameters of the internal 
(Ising) system J, p, and K are taken to be independent of 
temperature and bubble loop entropy is not explicitly in- 
cluded; the melting transition is therefore driven only by 
the difference in bending rigidity. The partition function 
z = E{ CTl } I Hi ^-e~ pH ^^ in transfer kernel form is 

TV 



{cr i= ±l} i=l 



(6) 



where \V) = (e^/ 2 , e~ p / 2 ) is the end vector, which enters 
to account for the free chain boundary conditions, and 
P(0^, f2j+i) is the transfer kernel given by 



P 



K J 4[cOS7i— l] + J + fl+K gKAB[cOS7i— 1] — J 

e ^AB [cos 7j — 1]— J g K B [cos 7j — 1]+ J— ix— K 



(7) 



with cos7i = tj • tj +1 . The A and B states form 
the canonical base, \A) = | + 1) = (1,0) and \B) = 
| — 1) = (0,1). Thanks to the rotational symmetry 
of the bending energy we can again integrate out the 
chain, leading to an effective Ising model with a free 
energy, iii^o, containing renormalized parameters: Z = 
e- {N - 1)r °E{a t} e-^.oM where H Ifl = H^K^p) 
with J a = J — [Gq(k a ) + Gq(kr) - 2Gq{k A r)}/4: and 
K Q = K-[G (ka)-G (k b )}/2 = K-AG$ B /2 the renor- 
malized Ising parameters, and To = [Gq(ka) + Gq(kr) + 
2Go(k j 4b)]/4. These parameters depend on chain rigidi- 
ties through Go(k) which is the free energy of a single 
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joint (two monomer) subsystem with rigidity k: Gq(k) = 
- M/ l§ exp[K(cos(0) - 1)]} = k - In[sinh(K)//c] . The 
renormalized quantity 2L = 2(/x + -Ko) corresponds to 
/3-Dr in the continuum model. If the bending free energy 
gain in opening one link, AGq B , is greater than the in- 
trinsic energy cost, 2(/j,+K), of opening an interior bond, 
then Lq becomes negative, signaling a change in stability 
of the A and B states. In the limit of high ka and kb, the 
entropic contribution dominates: Gq(k) sw —So(K)/kB 
ps ln(2«;). The discrete model then reduces to the con- 
tinuum one and AGq B w Iii(ka/ Kb) corresponds to the 
correction appearing in /3-Dr,. The difference between L$ 
and \i when K =^ creates an end-interior asymmetry 
that plays an important role in finite size effects. 

The Ising partition and correlation functions are 
obtained using transfer matrix techniques. The 
eigenvectors, |0, ±), and the eigenvalues, Ao,± = 
e Jo - r °[cosh(L ) ± (sinh 2 (L ) + e'^ ) 1 / 2 ], of the effec- 
tive Ising transfer matrix allow us to calculate the (di- 
mensionless) free energy per Ising spin of the coupled 
system, F = - InZ/N, where Z = £ T=± (V|0, t ) 2x o,t 
The average of the internal state variable is (c) = 
(J2i=i <fi/N) = —dF/dfi, from which the fractions of 
A and B links, <fA = (1 + (c))/2 and <pb = 1 — fA, can 
be derived. The melting temperature T m is then defined 
by (fiB{T m ) — 1/2. When N — > oo, (c) gets simplified to 



d In Ao,. 
9/i 



sinh(Lo) 



[sinh 2 (i ) + e~ 4J o 



1/2' 



(8) 



If Lq vanishes at a temperature, T™, sufficiently low for 
the cooperativity, or loop initiation, factor, a = e _4,/o , 
in the denominator to be small, then the system will un- 
dergo a melting transition: (c)oo will sharply cross-over 
from +1 for T < (pure A state) to -1 for T > T™ 
(pure B state). Contrary to previous Ising-type mod- 
els, the melting temperature is not put in by hand but 
emerges naturally from Lq = 0. In Eq. {8} e~ AJ ° de- 
termines the width of the transition region: ATf£ = 
2\d{c) OQ /dT\~l a (2M^) 2 /M)exp[-2J (T-)]. In 
the limit N,i — * oo , the influence of end-monomers dis- 
appears and the Ising correlation function reduces to 
((o- t -(c) 00 )(a l+r -(c) oc )) -> (l-(c)^)exp(-r/£i), where 
£i = — ln _1 (Ao,-/Ao.+) is the Ising correlation length, the 
typical size of B (A) domains below (above) T m . 

The tangent-tangent correlation function, (tjj_ ti+r), is 
obtained using the full transfer kernel method [ljj , which 
requires solving a spinor eigenvalue problem: P\^} = 
A|*), where |#(0)> = (*+i(0), We find eigen- 
values, A; T , labeled by Z = 0, . . . , oo and r = ± with the 
same form as for I = given above, but now where Go 
in the renormalized parameters is replaced by Gi(k) — 
k — m[ft'(d/ftd/c)'(sinh(fi;)/fv)] [2(j. The eigenspinors are 
|*j,m;r(0)) = V^r^m (ft) \l, t) with Y Zm (fi) the spherical 
harmonics. The transfer kernel can be expanded in terms 
of the eigenspinors P = J2i m T ^i,r|*/,m;r)(*i,m;r| and 



then be used to calculate the correlation functions. In the 
limit N,i — > oo, the expression for the tangent-tangent 
correlation function simplifies to 



(t< • t, 



+rl 



N. 



E 

T=± 



<l,r|0,+) 2 exphr/^. r ] (9) 



which reveals the importance of the two persistence 
lengths, Ci± = — In 1 (Ai : ±/ Aq,+) (units of £ ). It is 
not possible to extract from (tj • ti+ r ) one unique length 
for the whole range of T, because the weights associated 
with £f ± strongly vary with T. 

We now compare the discrete model predictions 
with experiment [l[ by focusing on the melting pro- 
file, Lps{T), of a synthetic nonalternating homopolynu- 
cleotide, polydA-polyrfT. A typical profile for a ho- 
mopolynucleotide has a sigmoidal shape characterized by 
T m and AT m . Of the five independent parameters that 
appear in the theory when K = 0, three are determined 
experimentally (polymerization index N and bending 
moduli assuming kab — ka [17]) and two (/x and J) are 
determined by fitting the model to experiment, hence 
T m = T m (p,,J;N,KA,K B )- Figure UK shows tp B {T) for 
DNA of molecular weight 1180 kDa [lj. From the known 
persistence lengths we obtain ka = 147 and kb = 5.54 at 
300 K. The solid line corresponds to our model fit with 
fb = 4.46 kJ/mol and J = 9.13 kJ/mol, leading to T m = 
326.4 K. The fitted value for 2/i is close to the experimen- 
tal energy of 10.5 kJ/mol needed to break an A-T link 0]. 
The value for J ~ 2jl is also consistent with the idea 
that destacking energy makes the dominant contribution 
to DNA stability The renormalized cooperativity pa- 
rameter at T m is Jo = 11.5kJ/mol > J. The model fit 
thus leads to parameter values in accord with experiment 
(in reality, the fitted values of fl and J implicitly com- 
pensate for effects like loop entropy explicitly left out of 
the model 

0,113). InFig.UJi, the curve corresponding to 
N — > oo is shown for the same parameter values. In this 
case, <pb(T) is given by Eq. JHJ) and is obtained from 
L a = 0: k B T™ ~ 2(fi + K)/ln(k A /KB)- In this limit the 
transition width is non-zero, due to the finite coopera- 
tivity parameter: AT~ oc ^\T™ ) = 2exp[-2J (T~)]. 
Since £i(T^) ~ 2000, typical helix and bubble domains 
are flexible within a small window about Tf£ . When 
N decreases, the width increases (Fig. QJi) roughly as 
AT m {N) - ATf£ ~ JV _1 . Hence even for a long poly- 
mer (N ~ 10 3 ), finite size effects are non-negligible, 



in agreement with experiments 211 ] . Then the nature 
of end monomers becomes important, as confirmed in 
Fig. [TK, which shows how short chains begin to melt by 
end-unwinding at lower temperatures (the trend that T m 
increases with decreasing N will likely be reversed when 
loop entropy is included @, [I3])- Our model predictions 
for experimentally accessible A-T_pair quantities are in 
agreement with accepted values 0, Q: i) a « 10~ 7 at 
T m ; and at physiological temperature ii) an interior sin- 
gle base-pair opening probability of 10~ 6 with a bub- 
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homopolynucleotide solutions 23[, in qualitative agree- 
ment with Fig. [TJd. Incorporating bending rigidity into 
DNA denaturation models thus allows us to make explicit 
predictions for both melting profiles and DNA mean-size 
dependent quantities. It will be of great interest to both 
probe such effects by carrying out experiments on DNA 
homopolymers and other biopolymers undergoing helix- 
coil transitions and extend our theory to heteropolymers, 
mechanical denaturation, and DNA dynamics. 
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FIG. 1: (a) Melting curves (Fraction of broken base- 
pairs vs. temperature) for polyrfA-polydT (data for 0.1 
SSC (= 0.015M NaCl +0.0015M sodium citrate, pH 7.0), 
N =1815 0]). The solid line represents the theoretical result 
for ft = 1.64 k B T m , J = 3.35 k B T m (T m = 326.4 K). The 
broken line corresponds to N — > co (same parameter values). 
Lower inset: melting curves for N = 100, 500, 1815, oo (in de- 
creasing order at low temperture, T < 326 K); Upper inset: 
Model results for the shift in transition width AT m — AT^ 3 
vs. polymer length [24| . (b) Temperature variation (N — ► oo) 
of the Ising correlation length, £i; persistence lengths of the 
coupled system, £^ ff and ± ; and of the pure chains, £ A B 
(in units of £q). At , the effective persistence length, 
rapidly crosses over from £ A to £, B . 



ble initiation barrier of 17fceT , and iii) a free energy 
of 0.l8kgT for breaking an additional base-pair in an 
already existing bubble. 

In contrast to purely Ising- type models, included in 
the predictions of our theory are mechanical and struc- 
tural features of the fluctuating chain, such as persis- 
tence length or mean-square-radius, R = (R 2 ) 1 / 2 with 
R = £q JN t i . From the expression for R in the limit 
N — > oo, we can define an effective persistence length, 
e eB : R 2 * (2£ 2 Q N)e eS = (2^A0£ T=± (0,+|l,r) 2 ^ r . 
Due to the coupling between bending and internal states, 
the respective weights (0, +|1,±) 2 associated with each 
correlation length change abruptly at T m (cf. Fig. [lb). 
Since the transition is very abrupt, it should also be pos- 
sible to observe it experimentally by measuring the radius 
of gyration by tethered particle motion [g^], light scat- 
tering, or viscosity experiments. For instance, since the 
relative viscosity is proportional to i? 3 , it should clearly 
exhibit an abrupt thermal transition. Such a transition 
has indeed been observed for the viscosity of synthetic 
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